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ABSTRACT 

Using the optical data from the Wide component of the CFHT Legacy Survey, and new 
ultraviolet data from GALEX, we study the colours and specific star formation rates 
(SSFR) of ^ 100 galaxy clusters at 0.16 < z < 0.36, over areas extending out to radii of 
r ~ 7 Mpc. We use a multicolour, statistical background subtraction method to study 
the galaxy population at this radius; thus our results pertain to those galaxies which 
constitute an excess over the average field density. We find that the average SSFR, and 
its distribution, of the star-forming galaxies (with SFR> 0.7 Mq/yt at z ^ 0.2 and 
SFR> 1.2 Mq/ji at z ~ 0.3) have no measurable dependence on the cluster-centric 
radius, and are consistent with the field values. However, the fraction of galaxies with 
SFR above these thresholds, and the fraction of optically blue galaxies, are lower for 
the overdense galaxy population in the cluster outskirts compared with the average 
field value, at all stellar masses > 10^'^ M0 and at all radii out to at least 7 Mpc. 
Most interestingly, the fraction of blue galaxies that are forming stars at a rate below 
our UV detection limit is much higher in all radial bins around our cluster sample, 
compared with the general field value. This is most noticeable for massive galaxies 
M* > lO^"'^ Mq; while almost all blue field galaxies of this mass have detectable star 
formation, this is true for less than 20% of the blue cluster galaxies, even at 7 Mpc from 
the cluster centre. Our results support a scenario where galaxies are pre-processed in 
locally overdense regions, in a way that reduces their SFR below our UV detection 
limit, but not to zero. 

Key words: galaxies: clusters: general, galaxies: evolution 



1 INTRODUCTION 

A lot has been learned about the star formation history of 
our Universe since early studies showed that the global star 
forma tion rate has declined by about a facto r of 10 since z ~ 
1 fe.g iLillv et"allll996l : [Madau et al.lll996l '). With decreas- 
ing redshift, the star formati on rate of star-forming galajc- 
ies h as been decreasing (e.g. iBell et al I I2OO5I : iNoeske et al.l 
120071). with an increase in the number of passive galax - 
ies (|Pozzetti et all I2OI0I '). The study by iBell et all (|2007l ') 
showed that a transformation from blue galaxies to red ones 
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is needed, so that the stellar mass of today's blue galaxies 
is not overproduced. 

The origin of this decline of the star formation rate and 
transformation remains unclear. One possible interpretation 
is that it is linked to the changing environment of galaxies. 
It is known that in dense regions such as the cores of galaxy 
clusters, the population is dominated by gala xies with red 
colours and low average star formation rates (|Balogh et al 



2004| : IWeinmann et al.l20 06': Haines ct al.'2006|- |Kimm et al 
20091 ). Since, under the hierarchical paradigm, galaxy clus- 
ters grow by accreting galaxies that are generally star- 
forming, a transformation must have happened to quench 
the star formation. However, how and where this happens 
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remains elusive. Studies have shown that the suppression 
of star formation is not re st ricted to clus t er cores. Fo r ex- 
amole. iGoinez et al.l (|2003l ). iBaloeh et aP (|l997l . Il998l l and 
iLewis et al.l (|2002 ) detected a lower fraction of star- forming 
galaxies relative to the field beyond ~ 2 virial radii. 

Despite the low star-forming fraction in clusters, the de- 
tection of the change of star formation rate within the star- 
forming population itself is still an unsettled issue. Some 
studies found that the specific star formation rate of star- 
forming gala xies is the same in the lowest and highest den- 
sity regions dPeng et al.l l2010l ) , in groups and in the field 
llMcGee et al.ll2011 ^ and in different regions in a z; ~ 0.2 su- 
pe rcluster (iBiviaiip ct al. 20l3); but (for example) the study 
bv IVulcani et aLn 20101 ) found a lower SFR of star- forming 
galaxies in 2 ~ 0.5 clusters than in the field. 

An interesting place to look for the transformations is 
the outskirt regions of clusters. Simulations have shown that 
clusters are surrounded by large-scale structures such as fil- 
aments and sheets , and galaxies are accreted mainly alori g 
these structures (iBond et al.l 1 19961 : IColberg et all Il999l ). 
These accretion zones are where galaxies reside before they 
reach the cluster cores, and thus might be an attractive 
candidate for where the transformation of galaxies hap- 
pens. However, studying the infall region of clusters is very 
difficult due to the low density contrast with the fore- 
ground/background field, and thus the exploration has only 
begun relatively recently. Studies focusing on two interme- 
diate redshift clusters and the Shapley supercluster found 
evidence of obscured star formation and the tr ansformation 
of spiral galaxies into SO i n infalling groups ( e.g. iGeach et al.l 
I2OO6I : iMoran et al.ll2007l : iHaines et a.l.ll201li ). 

In this paper, we use data from our CFHT Legacy Ul- 
traviolet Extension (CLUE), to study the star formation 
properties of a large sample of clusters, from the cluster 
core out to ~ 7 Mpc (the typical virial radius of clusters in 
our sample is ~ 1 — 2 Mpc). We examine these properties as 
a function of stellar mass, as it has become clear that stel- 
lar mass is one of the key parameters that d etermines the 
properties of galaxies (jKauffmann et aLll2003l ). In Section |2l 
we describe the data sets, and the cross-matching of the op- 
tical and UV catalogues. We describe our cluster and field 
samples, and the background subtraction procedure in Sec- 
tion [3] The stellar mass and SFR estimates are described in 
Sectional We present our results on the star formation rate 
and fraction of blue/red galaxies as a function of cluster- 
centric radius in Section (5] We discuss the implications of 
our results in Section |S] and conclude in Section [7] 

We assume a cosmology with Q.m ~ 0.3, JIa ~ 0.7 and 
_ffo=70 km Mpc~^ throughout. All magnitudes are in 
the AB system unless otherwise specified. 



2 DATA 

The data we used in this study are from the Wide component 
of the CFHTLS optical survey and our extended GALEX 
coverage over the CFHTLS fields. We describe them sepa- 
rately below. 



2.1 Optical Data 

The cluster sample we used in this study is detected from 
the CFHTLS Wide survey, which is a joint Canadian and 
French imaging survey in u* , g' , r' , i' and z' filters using the 
wide-field imager MegaCam. The survey is now complete, 
covering a total of 171 square degrees, composed of single 
pointings each with a field of view of 1x1 square degree. The 
total exposure time are 6000s in u* , 2500s in g' , 2000s in r', 
4300s in i' , and 7200s in z' . For extended sources, the 100 
per cent completeness limit is 25 mag in u* , ~ 25 mag in 
g' , ^ 24 mag in r', ~ 24 mag in i' , and ~ 23 mag in z' . More 
details of the data sets are provided in lLu et al.1 (I2009D . and 
thus here we only point out a few improvements. 

One of the improvements is that we used the stellar lo- 
cus of stars from SDSS to calibrate the colours of each indi- 
vidual CFHTLS pointing. The {g'-r') and (r'-z') colours of 
stars in each pointing are forced to agree with a reference set 
that is calibrated against S DSS stars but remain s on the na- 
tive MegaCam system (see lCilbank et aLll201ll for details). 
This reduces the scatter in colour when we stack a large 
number of clusters together. To calibrate the magnitude in 
each filter, we hold the r'— band magnitude fixed, and adjust 
the zeropoint of the other filters, on a field-by-field basis, 
to match the colour distribution of the stellar calibration 
set. This is because the r'— band magnitude from Terapix, 
calibrated internally with respect to some reference point- 
ings through overlap regions, has small variation between 
pointings, and is more accur ate than th e calib ration against 
2MASS, as described in Gil bank et al] (120111 ). The disper- 
sion of the final calibrated magnitude among all pointings is 
~ 0.01 mag, similar in all filters. The other improvement is 
that we excluded regions around bright stars, which would 
lead to photometry with larger uncertainty and an underes- 
timate of the galaxy number density in these regions. This 
second effect has the most impact when we study the prop- 
erties of galaxies far from cluster centres, where the density 
is low. 



2.2 GALEX Data 

In Cycle 5 of the GALEX Guest Investigator Programme, we 
proposed extended GALEX NUV coverage over the whole 
CFHTLS Wide fields (GI5-28fl to the same depth as the 
Medium Imaging Survey (~1500s). The data collection is 
only partially complete, currently covering approximately 
80 square degrees of the sky (about 50 per cent of the 
legacy field). In addition, we include existing archival data 
(57 pointings) over the CFHTLS Wide fields with exposure 
time greater than 1500s, which brings the total area with 
GALEX coverage to 110 square degrees. For multiple ob- 
servations taken at the same position, the one with deeper 
exposure is used. Figure [T] shows all the data we have in 
hand that are used in this work. 

We make use of the photometry pro vided in the NUV 
catal ogue measured using SExtractor llBertin fc ArnoiiU 
WQft) by the standard GALEX pipeline ( Morrissev et al. 
20071 ) . The standard pipeline uses a sophisticated method 



^ The full catalogues will be published online once the survey is 
complete; in the meantime the current catalogues are available 
upon request. 
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Figure 1. The current GALEX NUV coverage over the four 
CFHTLS Wide fields, including both our GI data and archival 
data. Points show objects brighter than NUV=23.0 mag. 




Figure 2. An example of the regions with residual elevated back- 
ground in GALEX NUV images, which causes problems for both 
object detection and photometry measurement. The size of this 
region is about 0.15 x 0.15 square degree. The unit on the colour- 
bar is counts per second per pixel. 



to correctly estimate the background in low-count regions. 
However, in some regions, there is residual elevated back- 
ground, likely due to reflection from bright stars. An exam- 
ple is shown in Figure (2] These elevated backgrounds cause 
problems in both object detection and photometry measure- 
ments. Therefore, we mask out these regions when perform- 
ing the analysis. The total area masked out is about 0.3 
square degrees, which is insignificant compared to the whole 
coverage; however it can become more important when con- 
sidering individual clusters. 



The fuU width half maximum (FWHM) of the GALEX 
point spread function (PSF) is ~ 4 — 5 arcsec, and the field 
of view is circular, with a radius of ~0.6 degree. To avoid the 
slightly degraded astrometry and photometry near the edge 
of each pointing, we only include objects within a radius 
of 0.58 degree from the field centre. As there are duplicate 
objects in the overlaps among adjacent pointings, we keep 
the one that is closest to the centre of the pointing it comes 
from. 

We use MAG_AUTO as the total magnitude. We correct 
for Galactic extinction using the relation between ex tinction 
and the reddening determined bv lWvder etHI ll2007h. where 
the reddening E(B-V) is calculated using the Schlegel et al-l 
(1998) dust map. Since the PSF of GALEX data is much 
broader that of the CFHTLS data, when measuring NUV- 
optical colours we use mag_AUTO in both catalogues, instead 
of aperture magnitudes as for the optical-optical colour mea- 
surements. We estimate the NUV magnitude uncertainty 
using duplicate detections from overlapping pointings (note 
this gives an upper limit of the uncertainty because the over- 
lapping regions are on the edge of each pointing, where the 
photometry is most uncertain). The standard deviation of 
the magnitude difference between duplicate objects is plot- 
ted in Figure [3] as a function of magnitude. To estimate 
the completeness of our GI5 data, we compare the number 
counts from our GI5 data with that from the Deep Imag- 
ing Survey (DIS) in the XMM-Newton Large-Scale Strcture 
(XMMLSS) Survey fields. In the top panel of Figure H the 
red dotted histogram shows the number count as a function 
of magnitude (normalized to per tile) from our GI5 data, 
and the black solid line shows the number count from the 
deeper data in the XMMLSS fields. In the bottom panel, we 
plot the ratio of the number counts from GI5 data to that 
from the XMMLSS fields, which shows that the ~ 80 per 
cent completeness of our GI5 data is about NUV=23.0 mag. 
Therefore, for our analysis here, we only consider objects 
brighter than this limit. 

2.3 Cross-matching Catalogues 

Since the CFHTLS data are much deeper than our GALEX 
data, essentially all NUV detected sources should have a 
match in the optical catalogue, except for extremely blue, 
optically faint galaxies {NUV — u* ^ —3 and u* ^ 26) that 
are not of interest here. Therefore, for each NUV source, we 
search the optical catalogue (without any magnitude cut) for 
matches within a radi us of 4 arcsec (the opt imal matching 
radius as discussed in iBudavari et al]l2009l ). In ~ 40 per 
cent of cases, there are multiple candidate optical matches 
for one NUV source. To deal with this, we take the closest 
match, unless there is a second candidate within 1 arcsec 
of it. This occurs for about 35 per cent of the cases with 
multiple candidate matches (so 14 per cent of all galaxies). 
In such cases, we use the colours to help identify the most 
likely counterpart, under the assumption that the most likely 
match is the one that has the most common colour for a 
galaxy of its magnitude. This procedure is described in the 
Appendix. 

We restrict the final matched catalogue to objects that 
are optically flagged as galaxies that are in the non-masked 
region (both around bright stars and regions with elevated 
background in the GALEX data, see Sections 12.11 and 12. 2p . 
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Figure 3. The standard deviation of the magnitude differences of 
dupUcate objects measured from overlapping regions in GALEX 
NUV images, as a function of NUV magnitude. 



with SExtractor flag<=3 (i.e. excluding objects close to im- 
age edge, with corrupted aperture or with at least one pixel 
saturated) in all filters {u* , g' ,r' ,i' , z' and NUV). 



3 SAMPLES AND BACKGROUND 
SUBTRACTION 

3.1 Cluster Sample 

Our clusters are detected from the optical CFHTLS data de- 
scribed above in Section [2m The overdensit y of red-sequence 
galax ies was used as tracers of clusters (|Gladders fc Yeg 
2OO0I) . The dete ction procedure was described in full detail in 
Lu et al.l (|2009l ). From the 171 square degrees, we detected 
~200 clusters with Nred,m'+2 > 12 in the redshift range 
0.16 < 2: < 0.36, where Nred,m'+2 is the number of red- 
sequence galaxies brighter than m* +2, within a radius of 0.5 
Mpc from the cluster centre. The richness Nred.m*+2 ~ 12 
roughly corresponds to a mass of 10^* Mq (see lLu et al.ll2009l 
for details). In Figure [5] we plot the richness vs. redshift of 
the clusters in our sample. Because the GALEX data are 
only available for about 65 per cent of the whole CFHTLS 
fields, for our analysis we use a subset of 112 clusters with 
GALEX coverage. To do the background subtraction as we 
will describe in the next Section, we divide our sample into 
two redshift bins, 0.16 < z < 0.27 and 0.28 < z < 0.36, and 
stack the clusters in each bin to increase the statistics. The 
number of clusters in the two bins are 43 and 69 respectively. 



3.2 Background Subtraction 

To study the properties of galaxies in clusters, it is important 
to separate them from the foreground/background galax- 
ies that are projected along the line-of-sight to the clusters 
of interest. We use the photometric redshift catalogue of 
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Figure 4. Top panel: number counts as a function of magnitude 
(normalized to per tile) from our GI5 data (red dotted histogram) 
and the deeper data in the XMMLSS fields (black solid line). 
Bottom panel: completeness of our GI5 data as a function of 
magnitude. In this work, we limit our sample to those brighter 
than NUV=23.0 mag, which is the ~ 80 per cent completeness 
limit. 



the C FHTLS Wide survey (|llbert et al.ll2006l : ICoupon et all 
l2009h to reduce the contamination by removing galaxies 
whose 2cr redshift uncertainties do not put them within 
Az = ±0.1 around the redshift of interest. We reduce the 
contamination further by using the photometric redshift cat- 
alogue to outline the region in the multi-colour space oc- 
cupied by galaxies at the redshift of interest, and exclude 
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0.25 0.3 

redshift 

Figure 5. Richness vs. redshift of the clusters in our sample, 
where the richness Af^ed,m*+2 is defined as the number of red- 
sequence galaxies brighter than m* + 2, within a radius of 0.5 
Mpc from the cluster centre. 



galaxies (including those without photometric redshift) that 
do not fall within the 99 per cent contour. 

Note that the redshift distribution of cluster members 
is much narrower than the typical uncertainty of the pho- 
tometric redshift. Therefore, a statistical background sub- 
traction is still needed, especially when we examine the out- 
skirt regions of clusters where the contrast with the field 
is low. Somewhat different from the usual background sub- 
traction method used in previous studies, we take advan- 
tage of our multi-band photometry and do the subtraction 
in multi-colour space, as we will describe in detail below. 
Our technique allows the removal of objects which do not 
have plausible colours to fit any typical SED at the redshift 
of interest. This will be important when we turn to models 
to convert observed quantities to physical ones. This multi- 
colour subtraction technique is somewhat analog ous to fit- 
ting empirical photometric redshifts (e.g., iBolzonella et al.l 
I2OOOI ). 

To construct the fore/background contamination sam- 
ple (referred to as the background sample from now on), we 
do exactly the same as with the cluster sample, but replace 
the centre of each cluster with a random position generated 
from the same patch (Wl,2,3,4) the cluster was detected 
fronfl. This way, for each cluster in the stack, there is a 
corresponding background sample with the same system- 
atics such as the imaging depth and overall number den- 
sity (which could vary from patch to patch due to different 
Galactic extinction and cosmic variance). 

Due to the masking (around bright stars and regions 



^ For the background sample, we mask out any known clusters, 
out to r = 3 Mpc, though this makes no appreciable difference. 



with elevated background as discussed in Sections 12.11 and 
I2.2|l and the gaps between GALEX tiles, the areal data cov- 
erage for each cluster is not 100 per cent. Therefore, we 
need to calculate, for each cluster, within the radius of in- 
terest, what fraction of the area is covered by valid data. 
Only clusters with a valid coverage greater than 80 per cent 
are included. We estimate the coverage fraction by using 
a high-resolution random catalogue, where points are ran- 
domly distributed over the footprint of our survey with an 
average density of 1 object per 10 square arcsec (a balance 
between computing efficiency and resolution requirement). 
We then weight each galaxy by this fraction, and by the 
number of clusters in the stack so that the number attached 
to each galaxy represents the number per cluster, i.e. 



/ / Ndata 



(1) 



where Ndata/ Nrand is the areal fraction covered by valid 
data and Ud is the number of clusters in the stack. 

With the two sets of samples in hand, one with cluster 
+ background counts, and one with just background counts, 
the background-subtracted net count of each galaxy, rii, is 
simply: 



■ Wcl + bg.i — 'Wbg,i, 



(2) 



where Wbg^i is the weight (as in Equation [1} of a back- 
ground galaxy, and Wci+bg,i is the weight of the nearest 
cluster-|-background counterpart to that background galaxy. 
The nearest counterpart is defined to be the galaxy in the 
cluster-f background sample that is closest to that back- 
ground galaxy in the space of NUV magnitude, {NUV — 
u*), {u* — g'), {g' — r'), (r' — i'), and {i' — z ) colours with 
equal weighl|f|. 

If the weight of the background galaxy, Wbg.i, is greater 
than that of its counterpart, Wci+bg,i, the resulting iii would 
be negative. In that case, we set Ui to zero, and subtract the 
"negative excess" from the next- nearest counterpart, and 
so on until the net excess becomes positive. In the end, we 
have a subset of galaxies from the cluster-|-background sam- 
ple that all have a positive number attached to them. These 
galaxies should essentially all be cluster members, because 
galaxies that have similar spectral shape as the background 
galaxies are effectively removed by doing the subtraction 
in multi-colour space. This way, not only is the total magni- 
tude distribution of cluster members recovered^, but colour- 
dependent quantities such as star formation rate can be de- 
rived. 

We perform the last two steps, constructing the back- 
ground sample and subtracting it from the cluster sample, 
100 times, and average the results over the 100 realizations. 
The error bars on the results presented throughout the paper 
are the standard deviation of the 100 realizations. Thus, they 
reflect the field-to-field variance in the background, which is 
the dominant source of statistical error in our measurements. 

As an example, in Figure |6] we show the result of one 



It makes no significant difference to the results whether or not 
we incorporate the colour uncertainties in this procedure. 
■* We have verified this by comparing the magnitude distribution 
of the cluster members obtained from our subtraction method 
and the more traditional subtraction method. 
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Figure 6. The spatial distribution of galaxies in the field of view (r<7Mpc) of one example cluster. Open cyan circles are all potential 
cluster members, as determined from their photometric redshifts. Solid red dots are galaxies that are left after the background subtraction, 
with NUV detected sources indicated in blue. The big black circle shows the 5 Mpc radius. 



realization of the background subtraction of one of the clus- 
ters in our sample. Open cyan circles represent all potential 
cluster members in the field of view before the subtraction, 
with photometric redshift pre-selection. Solid red dots are 
galaxies that are left after the subtraction, with UV detected 
galaxies indicated by the blue dots. The large black circle 
shows a radius of 5 Mpc from the cluster centre. The colour- 
magnitude diagrams (CMDs) of this cluster, in four radial 
bins, are shown in Figure [7] The symbols are the same as in 
Figure [S] with the additional black points indicating back- 



ground galaxies. The red-sequence is clearly visible in the 
core of the cluster. 



3.3 Phot-z Field Sample 

To compare the properties of galaxies in clusters with that in 
the general field, we construct a comparison field sample. We 
select, from the parent catalogue where clusters are detected, 
galaxies in the redshift range covered by our cluster sample 
using their phot-z. For the field sample, the phot-z of each 
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Figure 7. The colour-magnitude diagrams (CMDs) of the example cluster shown in Figure \E\ in four radial bins. The symbols are 
the same as in Fig ure [6l with the additional black points indicating background galaxies. The solid lines are the model CMD from 
iBruzual fc CharlotI |2003[). To guide the eye, the dashed line indicate colours 0.2 mag bluer than the red-sequence. 



galaxy is used in the calculation of the fc-correction (Section 
I4.1|l . In Figure [8] we show the number of galaxies in the field 
sample as a function of i'auto magnitude. Because of the large 
sample size, the statistical uncertainties on the field values 
in our analysis are negligible and thus are not plotted in the 
Figures shown in this paper. 



4 DERIVATION OF PHYSICAL QUANTITIES 
4.1 k-corrections 

After the background subtraction what we are left with are 
cluster members. But note, these members come from a 
number of clusters in a relatively wide redshift bin, stacked 
together. By doing so, the original redshift of the host clus- 
ter (determi ned from the re d-sequence galaxies as described 
in detail in ILu et all |2009| ) is no longer relevant. There- 
fore, for galaxies within each redshift bin, the median red- 
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Figure 8. Number of galaxies in the field sample, as a function of 
^'auto magnitude. Because of the large sample size, the statistical 
uncertainties on the field values in our analysis are negligible. 



shift of the bin is used when computing the fc-correction. 
The fc-co rrections are computed using kcorrect v4_1_4 
l|Blanton & Rowcis 20o3) ■ Because at the redshift of our 
cluster sample, the observed NUV magnitude probes the 
FUV wavelength, we reconstruct the rest-frame FUV mag- 
nitude from our six passband ph otometry using the public 
code of iBlanton fc RoweisI l|2007l ). The typical correction is 
about 0.6 mag. 



4.2 Stellar Mass 

We estimate stellar mass using the i-band mass-t o-light ra- 
tio, as a function of SDSS {g — r) colour, from iBell et ahl 
l|2003l ). We project our k-corrected rest-frame {g' — r') colour 
measured i n CFHTLS filters onto t he SDSS system, using 
the code o f iBlanton fc RoweisI (|2007h . Note that this conver- 
sion is provided for a diet Salpeter IMF, and therefore to be 
consistent with the Kroupa IMF used in the UV SFR esti- 
matio n, the resulting stellar mass is scaled by 0.7 (|Bell et al.l 
I2OO3I ) to con vert to the standard Salp eter IMF and then di- 
vided by 1.5 (|Brinchmann et al]|2004l ) to convert to Kroupa 
IMF. Because the mass-to-light ratio we used is only a func- 
tion of colour, we have a well defined stellar mass complete- 
ness as a function of colour. In Figure [9] we show for our 
higher redshift cluster sample, the colour vs. stellar mass. It 
shows that for red galaxies we are complete down to lO'' * 
M0 , as indicated by the vertical solid red line. We only con- 
sider galaxies above this stellar mass in the analysis. 



4.3 Star Formation Rate 

The UV luminosity we observe cannot be directly converted 
to a star formation rate because of the existence of dust 
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Figure 9. The colour vs. stellar mass plot for galaxies in our 
higher redshift cluster sample. The mass-to-light ratio is a func- 
tion of (g' — r'), thus for the red galaxies we are complete down 
to 10^'* Mq, as indicated by the vertical solid red line. 



around star-forming regions. The dust will absorb the UV 
emission and re-emit it at a longer wavelength. Therefore, 
it is important to correct for the effect of dust absorption, 
before we can infer the star formation rate from the UV 
luminosity. 

The total IR to UV luminosity ratio (Ltir/Luv) is 
a reliable est imator of t he dust extinction in star-forming 
galaxies (e.g. lBuat|[l993 ). We do not have IR information 
for our cluster sample. However, it has been shown that 
there is a correlation bet ween Ltir/Luv and the slope of 
the U V continuum (e.g. iMeurer et al.l 19991 : ICortese et al.l 
l2006h . We use this to estimate Ltir/Luv , and eventually 
the dust extinction. 

If fit by a power law, the UV continuum 
can be described as /(A) oc A'', where / is 
in unit of ergs cm~^ (|Calzetti et al.l 

1 1994) . Therefore, the slope /3 can be estimated as 
P = (log/2 — log/i) / (logA2 — logAi). In the redshift range 
of our cluster sample, the rest-frame wavelength range 
used to measure the UV slope, 1250A< A < 2600A, 
is shifted to the bandpasses of NUV and u* . There- 
fore, in our case, we reconstruct the rest-frame 
FUV and NUV from our observed multi-band pho- 



tometry again using the code of iBlanton fc Roweia 
(2007). Converting the AB magnitude to flux using 



m{AB) 



-2.51og/, 



.6 



-2.51og(/AA7c) - 48.6, 



where is in unit of ergs s ^ cm ^ Hz ^, we have: 
P = -0.4 (m^* - mjvt/v) / (logA2 - logAi) - 2 = 
2.2{mFuv — rriNuv) — 2, where mpuv and vtinuv 
are in the rest-frame. 

It has been pointed out that the relation between the 
UV slope and Ltir/Luv is different f or star-bu r sting galax- 
ies and normal star-forming galaxies (|Bellll2002l : rKong et al.l 
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Figure 10. SFR vs. apparent NUV magnitude, at the two redshifts respectively. The insets show the cumulative SFR distribution of 
galaxies ~0.2 magnitude brighter than the apparent NUV magnitude limit (NUV=23 mag). Different histograms are for different stellar 
masses, and the vertical lines indicate the adopted SFRuv limits for our analysis. Depending on the stellar mass, ~ 60 — 90 per cent of 
the galaxies near the NUV detection limit is below the SFR limits we adopt (0.7 MQ/yr at z ~ 0.2 and 1.2 MQ/yr at z ~ 0.3). 



l2004l : ICortese et aLlbOOd ) . Since most of the galaxies we are 
after here are nor mal star-forminR gala xies, we adopt the 
relation derived by ICortese et all ()2006l )fea. 5) for normal 
star-forming galaxies, which is: 

log{LTiR/LFuv) = (0.7 ± 0.06)/3 + (1.3 ± 0.06). (3) 

In the study bv lCortese et al.] (|2008l '). the relation be- 
tween Ltir/ Lpuv and dust extinction A(FUV) is derived 
using SED fitting for galaxies with different ages. There- 
fore, for galaxies off the red-sequence (defined as those with 
[g' — r') and (r' — i') colour bluer by 0.2 magnitude than 
the red-sequence), we adopt the relation for young star- 
forming galaxies (their r > 7 Gyr model); and for galaxies 
on the red-sequence we adopt the one for an older popu- 
lation (their r ^ 5.4 Gyr mode l) (see their table 1). As 
suggested by lCortese et al.l (|200^ . for galaxies redder than 
{FUV — NUV) = 0.9, their extinction is assumed to be the 
same as galaxies with {FUV - NUV) = 0.9. 

The dust-corrected magnitude is turned into a lu- 
minosity, and the st a r form ation rate is then estimated 
using the iKennicuttI ll 19981) relation. The factor of 1.5 
Brinchmann et alT2004 ) is to c onvert from S alpeter IMF 



Salpeteij|l955l ) to Kroupa IMF (|Kroupall200ll ) 



SFRvv{MQ/yT) = 1.4xlO"^*L,(ergss"' Hz"') x 10'' '''*Vl-5- 

(4) 

Because of the above described conversion between 
NUV magnitude and SFR, there is a scatter between the 
observed NUV magnitude and the derived SFR, as shown in 
Figure [TO] (small dots). Thus, the limiting NUV magnitude 
of our sample does not correspond to a precisely defined SFR 
limit. We can estimate the SFR limit by looking at the SFR 
of galaxies near the limiting apparent NUV magnitude. The 
insets in Figure [lO] show the cumulative SFR distribution of 
galaxies ^0.2 magnitude brighter than the NUV magnitude 



limit (NUV=23 mag) at the two redshifts respectively. Dif- 
ferent histograms are for different stellar masses. At z ~ 0.2, 
depending on the stellar mass, ~ 60 — 90 per cent of these 
galaxies near the NUV detection limit have SFR below 0.7 
M@/yr. Therefore, we adopt 0.7 M^/yr as the SFR detec- 
tion limit for our NUV magnitude limited sample at this 
redshift. Similarly, at z ~ 0.3, we adopt a SFR limit of 1.2 
M0/yr. These limits are indicated as the vertical lines in 
the insets. The use of a different limiting SFR at each red- 
shift bin means we must be cautious when interpreting trend 
with redshift as evolution. Our main purpose, however, is to 
compare the field and cluster populations at each redshift. 



5 RESULTS 

In this section, we present our findings on the specific 
star formation rate (SSFR), fraction of star-forming and 
blue/red galaxies of the cluster sample as a function of 
cluster-centric radius, and compare them with that of the 
field sample. Note that, as discussed before, we divide our 
sample into two redshift bins, and thus the limiting SFR will 
be different for each bin. 



5.1 SSFR vs. Stellar Mass 

It has been suggested that star-forming galaxies form 
a tigh t sequence in the SSFR vs. stella r mass (M*" 



plane (iNoeske et all l2007l : ISalim et all |2007| ) (but see e.g. 
ICowie fc Bargeijl2008l where a much larger scatter is found) , 
where SSFR is defined as SFR/M*. In Figure [TT] we show 
one realization of the SSFR vs. M» of our stacked clusters at 
two redshifts, in the top and bottom panels respectively. The 
black squares are blue cluster galaxies, while red squares rep- 
resent those that are on the red-sequence. The green crosses 
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Figure 11. One realization of the SSFR vs. stellar mass relation of our stacked clusters at two redshifts (top and bottom panels), for 
the two regions with the most contrasting density, the core (0 < r < IMpc) (left panels) and the outskirt (5 < r < 7Mpc) (right panels). 
The black squares are blue cluster galaxies, and the red squares represent those that are on the red-sequence. Galaxies in the phot-z field 
sample are plotted as green crosses. The slanted dashed lines indicate the UV detection limits. The short solid lines indicate the three 
stellar mass bins used in the analysis. 



are phot-z field galaxies. In the left and right panels, we 
present the relation in two regions with the most contrast- 
ing density: the core (0 < r < IMpc) and the outskirts 
(5 < r < 7Mpc). Note that our SFRuv detection limits (as 
discussed in Section |4.3|) are 0.7 Mq yr~^ at 2 ~ 0.2 and 
1.2 Mq yr~^ at z ~ 0.3 (as indicated by the slanted line). 
Therefore, the slope of the SSFR vs. M* of our sample is 
mostly driven by the detection limit, and does not refiect an 
intrinsic correlation between these quantities. For this rea- 
son, we do not attempt to fit the relation for our sample. 
Nonetheless, comparing the left panels with the right ones, 
it seems that the relation in the cluster cores and the out- 
skirts are both consistent with that in the field, showing no 
dependence on the environment. 

The ensemble average SSFR (total SFR/total stellar 



mass) of galaxies with SFR above our UV detection limit 
is shown in Figure [T^ at different cluster-centric radii, split 
into three stellar mass bins. Black squares, red crosses, and 
blue triangles represent the 9.8 < log]^o(^*/^0) < 10.3, 
10.3 < logio(M,/M0) < 10.7, and 10.7 < logio(M./M0) < 
11.2 mass bins respectively. For all stellar masses, the aver- 
age SSFR of the star-forming galaxies (above our adopted 
SFRuv limits) is approximately constant from the cluster 
core out to the outskirt region. The Spearman's rank cor- 
relation tests show that the observed relation between the 
average SSFR and the radius deviates from no-correlation 
hypothesis on 1 — 1.5 a level for the various stellar masses. 
Note that, our errorbars prevent the accurate determination 
of the rankings of the variables; if ranking the average SSFR 
randomly, the level of deviation varies from ^ 0.3 — 1.7 a. 
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Figure 12. The ensemble average SSFR (total SFR/total stellar mass) of galaxies with SFR above our UV detection limit, at different 
cluster-centric radii, split into three stellar mass bins. Black squares, red crosses, and blue triangles represent the 9.8 < logj^g(M,/M0) < 
10.3, 10.3 < logjQ(M,/MQ) < 10.7, and 10.7 < logjQ(Mt/M0) < 11.2 mass bins, respectively. The corresponding field values are 
indicated by the black solid, red dotted, and blue dashed lines respectively. The average SFR of the star-forming galaxies (above our 
SFRuv limits) is roughly constant from the cluster core out to the outskirt regions, and is consistent with the field value within the 
uncertainties. 




Figure 13. The distribution of the SSFR of galaxies with 10.7 < logio(M*/M0) < 11.2 at 0.16 < z < 0.27, normalized by the total 
stellar mass, at different radii. The green circles, black solid and red dotted histograms, and cyan squares represent the four cluster-centric 
radii, and magenta triangles represent the field values. For clarity, we split it into two panels with the field values plotted in both as 
references. Within the relatively large error bars, there docs not seem to be a significant difference between the distribution at different 
cluster-centric radii and the field. 
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It is also interesting to compare this with the field val- 
ues, which are indicated by the black solid, red dotted, and 
blue dashed lines for the three stellar masses respectively. As 
mentioned before the statistical uncertainties on the field 
values are entirely negligible, given the large sample size. 
Within the uncertainties, there does not appear to be a sig- 
nificant difference between the average SSFR of star-forming 
galaxies in clusters and in the field. 

We now further examine the distribution of the SSFR 
(instead of just the average) in each stellar mass bin. For 
clarity, we only show the one for the most massive galaxies 
(10.7 < logio(M,/M0) < 11.2) in the lower redshift bin 
(0.16 < 2 < 0.26), as the behavior of galaxies at other stellar 
masses and at the higher redshift is qualitatively the same. 
In Figure 1131 the green circles, black solid and red dotted 
histograms, cyan squares, and magenta triangles represent 
respectively the distribution (normalized by the total stellar 
mass) of the SSFR at four cluster-centric radii and in the 
field. For clarity, we split them into two panels with the 
field values (magenta triangles) plotted in both panels as 
references. The distributions appear to be consistent within 
our relatively large error bars. This implies that, above our 
SFRuv detection limits, either the quenching timescale is 
short so that the population of transition galaxies eludes 
detection, or that only a small fraction of the population 
above the limits are going through the transition and the 
majority of the star-forming galaxies are left alone. 

There re s ults ar e i n agreement w ith Peng et al. ll201Cf): 
iMcGee et al] l|201lD : iBiviano et all (|201ll ): I Wetzel et all 
l|201ll) in general, but in apparent conflict w ith the low SSFR 
found in clusters by IVulcani et al] (|2010l '). We will discuss 
this further in Section [S] 



5.2 Star-forming Fraction 

We now examine the fraction of star-forming galaxies out 
of the whole population as a function of environment. To 
obtain the total number of galaxies, we follow the same pro- 
cedure as for the optical-UV matched catalogue, except that 
we now include galaxies that are not NUV detected. 

In Figure [T4l the fraction of star-forming galaxies with 
SFR above our UV detection limit is plotted as a function 
of cluster-centric radius in three stellar mass bins, for the 
lower (left panel) and higher redshift samples (right panel). 
Symbols are for cluster galaxies at three stellar masses (as 
labeled in Figure llSp , and the lines are corresponding field 
values. For all stellar masses at both redshifts, the fraction 
of star-forming galaxies in clusters is much lower than that 
in the field, and is independent of the radius within the 
uncertainties. The very low star-forming fraction in clus- 
ters is partially due to our relatively high SFRuv detection 
limits. Note that, the stellar mass dependence of the star- 
forming fraction shown here (for both clusters and the field) 
may seem counterintuitive, but keep in mind that our star- 
forming fraction is calculated with fixed cuts in SFR instead 
of SSFR. Therefore, although the star-forming fraction is 
higher for low mass galaxies, their typical SFR is lower and 
thus more of them would be below our SFR cut. In any case, 
our focus here is the comparison between clusters and field 
in each redshift bin, not the stellar mass dependence. 

In the next section, we examine the behavior of the blue 



population that includes galaxies with SFR lower than our 
UV detection limit. 



5.3 Blue/Red Fraction 

To estimate the fraction of galaxies that are on the red- 
sequence, we refine the model slope of the red-sequence from 
our data directly. We stack clusters in each redshift bin 
to a central redshift (pa ssively evolved based an old SSP 
iBruzual fc CharlotI |2003| model) , and fit the red-sequence 
and the scatter, cr, around the fitted colour-magnitude rela- 
tion. 

We consider galaxies that are redder than 3o" below the 
fitted red-sequence as red galaxie^. The fraction of blue 
galaxies as a function of distance from cluster centres is plot- 
ted in Figure [15] for the two samples at 2 ~ 0.2 (left panel) 
and 2 ~ 0.3 (right panel), again split into three stellar mass 
bins (black squares, red crosses, and blue triangles). At both 
redshifts, the blue fraction decreases with stellar mass, at all 
radii. At fixed stellar mass, it is quite clear that at both red- 
shifts, the change of the blue fraction happens within ~ 3 
Mpc from the cluster centres, with no significant further 
changes beyond ~ 3 Mpc. Furthermore, the difference be- 
tween the blue fraction in the outer and inner most regions 
is smaller for the most massive galaxies compared to those 
of lower masses. 

It is interesting also to compare this with the field val- 
ues. We calculate the blue fraction for the phot-z field sample 
in a similar fashion as for the cluster sample, i.e. red galax- 
ies are defined as those that are redder than 3ct below the 
fitted field red-sequence. The resulting field blue fractions 
are plotted in Figure [15] as the black solid, red dotted, and 
blue dashed lines for the three stellar masses respectively. 
Not surprisingly, in the central regions of the clusters (r < 1 
Mpc) the blue fraction is lower than that in the field. The 
more interesting thing is what happens in the outer regions. 
Beyond r ~ 3 Mpc, the cluster blue fraction is generally 
lower than that in the field, perhaps surprising given the 
low density contrast at these radii. Only the most massive 
galaxies, at 2 ~ 0.2, reach blue fractions comparable to the 
field at these radii. Keep in mind that due to the nature of 
the background subtraction, what we detect is the "excess" 
over the field population (at r ^ 7 Mpc, the excess is on the 
~10 per cent level with respect to the density in the field). 
We discuss the implication of these results in combination 
with the other results presented above in the next Section. 



6 DISCUSSION 

Our results show that the fraction of blue, actively star- 
forming galaxies is significantly lower in dense environments, 
even far from cluster cores. However, no significant dif- 
ference in the distribution of SFR for the active popula- 
tion (with SFR above our SFRuv detection limits) is de- 
tected. Both the average SFR and the shape of the distri- 
bution are the same in the field, the dense cluster cores, 



^ We have also confirmed that if defining red galaxies by mi rror- 
ing the redder half of the fitted red-sequence [lu et alj|2009h our 
conclusions do not change. 
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Figure 14. The fraction of star-forming galaxies witii SFR above our UV detection limit, as a function of cluster-centric radius in three 
stellar mass bins. Blue triangles, red crosses, and black squares represent three stellar mass bins (in decreasing order, as labeled in Figure 
I15II , and the corresponding field values are indicated by blue dashed, red dotted, and black solid lines. The lower and higher redshift 
samples are shown in the left and right panels respectively. 
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Figure 15. The fraction of blue galaxies as a function of distance from cluster centres, for the two samples at z ^ 0.2 (left panel) and 
z 0.3 (right panel), again split into three stellar mass bins. Symbols represent cluster values, while lines represent field values for the 
corresponding stellar mass. The blue fraction decreases with stellar mass. At fixed stellar mass, the change of the blue fraction happens 
within ~ 3Mpc from the cluster centres, with no significant further changes beyond ~ 3 Mpc. See text for a detailed discussion. 



and overdense regions in the distant outskirts. This is in 
good agreement with a growing body of independent re- 
sults. Studying the local universe, several authors have 
shown that the distribution of colour and/or SFR for star- 
forming galaxies in the SPSS is independent of environmen t 
jBalogh et all 12004 iPeng et all I2OI0I : IWetzel et al] l201lf) . 



Recently, iBiviano et all (|201ll ) came to the same conclusion 
using 24/xm-derived SFRs around a. z — 0.23 supercluster; 
while the fraction of star-forming galaxy depends sensitively 
on environment, in a complex way, the correlation between 
SFR and stellar mass rem ains unchanged. At even higher 
redshifts of 0.3 < z < 0.5, iMcGee et atl l|201ll 'l also found 
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Figure 16. Fraction of blue galaxies that have SFR above our UV detection limit. The two panels are for two redshifts. Different symbols 
represent cluster galaxies with different stellar masses, and the lines represent field galaxies at corresponding stellar masses. 



the mean SFR of star-forming galaxies in groups to be inde- 
pende nt of environment. Notab ly, the results of iPeng et al.l 
(|2010l ) and iMcGee et all (j201ll ;) extend to SFR lower than 
we probe in this paper, and yet they still find no population 
of low-SFR galaxies in clusters. 

This is i n ap parent contrast with the work by 
IVulcani et al.1 (|2010l '). who measured SFR from [Oil] and 
24/im data in the ESO Distant Cluster Survey, reaching a 
SFR depth comparable to ours. They found that the aver- 
age SFR of star-forming galaxies in these z ~ 0.5 clusters 
is a factor of 1.5 lower than that in the field. One possibile 
explanation for the discrepancy with our data is that if the 
low SFR popula tion in clusters are sufficiently dusty (e.g. 
I Wolf et al.ll2009l ). they might not make it into our NUV- 
selected sample. Galaxies at the limiting apparent NUV 
magnitude in our sample and with extinction greater than 
~ 1 — 2 mag would not have made it into our final SFR 
limited sample. 

We can further gain insight into the population with 
SFR below our UV detection limit (possibly including the 
dusty ones) by comparing the fraction of star forming galax- 
ies among the optically blue population. In Figure[Tn]we plot 
the fraction of blue galaxies that have SFR above our UV 
detection limit. The two panels are for two redshifts. Differ- 
ent symbols represent cluster galaxies with different stellar 
masses, and the lines represent field galaxies at correspond- 
ing stellar masses. As expected, not all blue galaxies are star 
forming (at a rate higher than our relatively high UV de- 
tection limit). However, what is especially remarkable is the 
difference between the field and cluster galaxies. Compared 
to the blue galaxies in the field, a lower fraction of the blue 
galaxies in the clusters is forming stars at a rate above our 
UV detection limit. This indicates a larger population of 
blue galaxies with low SFR (lower than our UV detection 
limit) in clusters than in the field. Although, it is possi- 



ble that, as mentioned above, these blue g alaxies in clusters 
are dustier than t heir field counterparts (|Wolf et al.l l200sj : 
iHaines et al.|[201ll ). and thus would still be optically blue 
but suppressed in the UV. If indeed their SFR is reduced 
and yet they are still blue, it would require a quenching 
mechanism that does not completely shut down the SFR on 
short timescales. Further, this mechanism must affect galax- 
ies in a way that is independent of proximity to the cluster 
core, as this fraction is roughly constant with radius within 
the uncertainties. 

It is also notable that, at the outskirts (~ 7Mpc), de- 
spite the low average density contrast with the field, cluster 
galaxies have a lower blue fraction than that in the field (ex- 
cept, perhaps, for the most massive galaxies at z~ 0.2). This 
indicates that those galaxies that just arrived at the clus- 
ter environment already have their star formation quenched, 
supporting a "pre-processing" scenario suggested by severa l 
studies (e.g. Zabludoff & Mulchaev 1998: Balogh et al.ll2000l : 
.Balogh fc McGee. 2010: McGee et al.,.20 09'). The recent work 
by IWetzel et al.l l|201ll ) also found enhanced fraction of 
quenched galaxies out to ~ 10r2oo in SDSS groups/clusters 
and attributed it to satellite galaxies that have possibly been 
pre-processed in groups. To gain more i nsight into t his is - 
sue, we use the semi-analytic models of iFont et al.l (|2008l ) 
to construct the halo mass function of galaxies at different 
cluster-centric radii an d compare it with that of the field 
galax;ies. The .Font et al.. (2008) model is a recent version o f 
the GALFORM model jCole et al.ll2"oO0l : iBower et al.ll2006l ) 
which has been modified to include a more realistic treat- 
ment of environmental effects. Although this method still 
has some proble ms in reproducing the distribution of galaxy 
colours (e.g. iBalogh et al. 2009). the stellar mass functions 
agree reasonably well with observations to z ~ 5. To mimic 
our observations of the large scale environment around clus- 
ters, we construct a stacked model cluster by selecting all 
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Figure 17. Fraction of galaxies (9. 8 < log(M»,/M(7)) < 10.3) as a function of cluster-centric radius which reside in halos above a given 
mass in the semi-analytic models of lFont et al.l ||2008|) . Solid curves, and dashed lines indicate cluster and field values respectively. 



galaxies within ± 4000 km/s of the central galaxy of clus- 
ters {Mh > lO^**-^ Mq) in the Font model. The advantage of 
using the models is that the host halo mass of each galaxy 
is known. In Figure 1171 we show the resulting fraction of 
galaxies with stellar masses of 9.8 < log(M«/MQ) < 10.3 
as a function of cluster-centric radius which reside in ha- 
los above a given halo mass. The lower and higher redshift 
bins are presented in the left and right panels respectively. 
Solid curves present cluster values, and the dashed lines in- 
dicate field values (see legend for details). It shows that at 
around ~ 6 Mpc, the outermost region we probed with our 
data here, the fraction of galaxies residing in halos above 
10^^' Mq is slightly higher than that in the field. This could 
be part of the explanation of our results: the higher red frac- 
tion in the outskirts is due to the fact that more galaxies are 
in halos above group scale halos in the outskirts than in the 
field, provided that the quenching already started in those 
halos. However, it is not clear whether those galaxies resid- 
ing in groups in the outskirts are further processed in the 
proximity of the large-scale cluster environment. Further- 
more, because of the background subtraction, what we ob- 
serve is the properties of the cluster population in "excess" 
of the field. Thus, it is hard from our data to determine the 
exact nature of the excess we detected at large radii. A com- 
parison between groups in the outskirts and isolated groups 
in the field will provide more insight. Therefore, we defer 
a closer examination of groups in the outskirt regions to a 
future paper, where we use spectroscopic data to study two 
contrasting clusters. 



7 CONCLUSIONS 

In this work, we examined the star formation properties of 
a large sample of ~ 100 galaxy clusters at 0.16 < z < 0.36, 
from their cores out to ~ 7 Mpc, using the CFHTLS optical 



data and GALEX UV data. Our main findings are summa- 
rized below. 

(i) We found that the average SSFR and the distribu- 
tion of SSFR (of galaxies with SFR> 0.7 Me/yr at z ~ 0.2 
and SFR> 1.2 M0/yr at z ~ 0.3) show no strong depen- 
dence on the distance from the cluster centre within the 
error bars, and are similar to that in the field as well. 

(ii) The fraction of star-forming galaxies (with SFR> 
0.7 Me/yr at z ^ 0.2 and SFR> 1.2 Mg/yr at z ~ 0.3) is 
much lower in clusters than in the field. For cluster galaxies, 
this star-forming fraction is radius and stellar mass indepen- 
dent within the uncertainties, partially due to our relatively 
high SFRuv limits. 

(iii) Among the optically blue population in clusters, 
a lower fraction is forming stars at a rate higher than 
our SFRuv limits compared to the blue population in the 
field. This difference is larger for high mass galaxies, and is 
roughly independent of radius. 

(iv) The fraction of galaxies with blue colours is con- 
stant from ~ 3 Mpc out to ~ 7 Mpc; however, within 3 
Mpc, there is an abrupt decrease in this fraction towards 
the cluster core. This is present at both redshifts, at all stel- 
lar masses examined here, but more so for the least massive 
galaxies (9.8 < logi(,(M,/MQ) < 10.3). 

(v) Despite the low average density contrast with the 
field in the outermost region (r 7 Mpc), the blue fraction 
is lower than that in the field; with the exception of the most 
massive galaxies at z ~ 0.2. 

Our results imply that the excess population over the 
field in the outskirts of clusters is pre-processed, and at all 
radii throughout the clusters there is a population of blue 
galaxies that have their SFR reduced to below our UV de- 
tection limit but not to zero. This requires a mechanism that 
does not shut off the star formation completely, and works in 
a way that is independent of radius. Limited by our SFRuv 
detection limit, we cannot probe the distribution of the SFR 
of that partially-quenched population. With deeper data to 
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detect these galaxies and measure the change in the shape 
of the SFR distribution, it may be possible to put stronger 
constraints on the timescale on which the quenching mech- 
anism operates. 
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APPENDIX A: CROSS-MATCHING 
CATALOGUES 

As discussed in Section [2.31 in cases where there are one or 
more candidate optical matches within 1 arcsec of the closest 
match to a NUV source, we use colours to help identify the 
most likely optical counterpart. We do this using galaxies at 
fixed NUV magnitude. For a galaxy at a certain magnitude. 
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Figure Al. One example of how we use colour information to 
cross-match NUV sources with optical sources. Small dots are 
galaxies at 22.5 < NUV < 23.0 and have only one possible op- 
tical counterpart. There is an obvious concentration of dots in a 
specific region in this colour-colour space, indicating the colours 
the real matches are mostly likely to have. The four squares are 
four possible optical matches within 4 arcsec from a NUV source 
with NUV=22.5, numbered in the order of distance from the NUV 
source, with 1 being the closest one. It shows that the one that 
is the closest match in this case is less likely to be the real match 
than the second closest match, because it is located in a less dense 
region. 

the probability of it having a certain colour is not random. 
In Figure [KT\ the (u* — r') colours of the objects that only 
have one possible optical counterpart (small dots) in the 
magnitude range 22.5 < NUV < 23.0 are plotted against 
their {NUV — it*) colours. As we can see, the density of the 
dots in this colour-colour plane is not uniform; instead, most 
galaxies have colours of (it* — r') ~ 1 and (NUV — u*) ^ 
0, and we interpret this as the most probable colour for a 
galaxy of this magnitude. Therefore, when there are multiple 
possible optical matches for a NUV source, we take the one 
that resides in the region with the highest density as the 
real match. One example is shown in Figure lAll The four 
squares are four possible optical matches to a NUV source 
with NUV=22.5 mag, numbered in the order of distance 
from the NUV source, with 1 being the closest. The one that 
is the closest match is residing in a region where the density 
of the points is less than that of the second closest match, 
and thus in this case we take the spatially second closest 
match as the real match. In about 15 percent of the cases, 
the spatially closest match is different from the colour-based 
match. 
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